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J^ , Abstract. We present observations of continuum and molecular lines towards the protostar candidates IRAS 



O 



21307+5049 and IRAS 22172+5549. Single-dish maps in the ^^CO (2-1), C^**0 (2-1), Hi^CO+ (1-0) lines and 



850 fim continuum are compared with interferometric maps in the ^^CO (1-0) line and 3 mm continuum, and 

r>^ , with mid- and near-infrared images. Also, single-pointing spectra of the CH3C2H (6-5), (8-7) and (13-12) lines 

Q ■ observed towards IRAS 21307+5049, and H2 and [Fell] line emission observed towards IRAS 22172+5549, are 

^H ^ presented. A plausible interpretation of our data based on the continuum maps and spectral energy distributions is 

f^ . that two components are present: a compact molecular core, responsible for the continuum emission at wavelengths 

longer than ~ 25 /im, and a cluster of stars located close to the center of the core, but not spatially coincident 

with it, responsible for the emission at shorter wavelengths. The core is approximately located at the center of 

the associated molecular outflow, detected for both sources in the ^'^CO (1-0) and (2-1) lines. The cores have 

O ' masses of ~ 50 Mq, and luminosities of ~ 10'' Lq. The outflows parameters are consistent with those typically 

found in high-mass young stellar objects. Our results support the hypothesis that in these sources the luminosity 

^ ' is dominated by accretion rather than by nuclear burning. We conclude that the sources embedded inside the 

' cores are likely protostars with mass ~ 5 — 8 Mq . 
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1. Introduction tion of disks and coUimated bipolar outflows driven by 

the forming star, whereas in the coalescence scenario the 

Molecular outflows are a well known phenomenon associ- outflows must be poorly coUimated: it is thus important 

ated with the star formation process. Their origin is re- ^^ investigate at high angular resolution the presence of 

lated to the formation of proto-stellar disks, and they play coUimated outflows originating from massive protostellar 

an important role in the accretion phase of a protostar. In candidates in order to discriminate between coalescence 

fact, they increase the turbulence of the surrounding ma- g^-^^ accretion 
terial, hence contributing to support the molecular cloud 

against gravitational contraction and perhaps halting the Although many studies have been carried out on out- 
collapse onto the forming star. The role of outflows is espe- ^^^^ origmatmg from low-mass objects, only a hmited 
daily relevant in the context of the high-mass (M > 8 Mq) ^^^^cr of observations are avaflable for those associated 
star formation process. While low-mass stars are believed ^i^^^ high-mass young stellar objects (YSOs). Moreover, 
to form through accretion, the mechanism for the forma- ^o^* ^^ ^^ese have beenperformed at low angular resolu- 
tion of high-mass stars is still unclear. Two scenarios have ^ion (e.g. Zhang et al. IMD Beuther et al. [20023. 
been proposed: accretion (Behrend & Maeder |2U(J1|I . or In this paper, we investigate the structure of two can- 
coalescence of several lower mass protostars (Bonnell et didate massive protostars: IRAS 21307+5049 (hereafter 
al. [TM5l Stabler et al.EOOOIl- Theoretical models predict 121307) and IRAS 22172+5549 (hereafter 122172). These 
that the accretion process naturally leads to the forma- sources, located at a kinematic distance of 3.2 and 2.4 
kpc respectively (Molinari et al. .200211 . are part of an mi- 
Send offprint requests to: F. Fontani, e-mail: tial large sample of luminous IRAS sources. In a series of 
fontaiii@arcetri.astro.it Studies, Molinari et al. I|1996lll998alll998b[l2000|l selected 
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from the initial sample of 260 sources, a small number of 
likely massive YSOs thought to be in an early evolution- 
ary stage preceding the formation of a compact Hii region. 
Then, more recently, Mohnari et al. f2002') observed 11 of 
these sources at high angular resolution with the Owens 
Valley Radio Observatory (OVRO) in the HCO+ (1-0) 
line and in the millimeter continuum, and in the 3.6 cm 
radio continuum with the Very Large Array (VLA). For 
several of these sources, a compact molecular core was de- 
tected, and in several cases this core does not show emis- 
sion at 3.6 cm. Moreover, in a few cases, non gaussian 
wings were identified in the HCO"*" (1~0) line, suggesting 
the presence of outflows. These results indicate that in the 
OVRO sample there is a number of sources which present 
observational features expected for protostars. This is con- 
firmed by Fontani et al. (2004) who observed one of these 
sources, IRAS 23385+6053, with single dish telescopes and 
interferometers at various wavelengths: the authors have 
demonstrated the presence of a massive dusty core, which 
likely hosts a young stellar object in the protostellar phase. 

Two other sources of the OVRO sample, 121307 and 
122172, show the presence of a compact dusty core with- 
out centimeter continuum emission, with prominent wings 
in the HCO+ (1-0) line, and characterised by IRAS lumi- 
nosities of a few 10'^ Lq (Molinari et al. 2002). Hence, 
they represent likely massive YSOs, in an evolutionary 
stage preceding the formation of a compact Hll region. 
Recently, Zhang et al. I|2001|l and Mohnari et al. (2002) 
have detected towards both sources a compact outflow in 
the I^CO (2-1), 13C0 (2-1) and HCO+ (1-0) lines. Here 
we present the results obtained from single-dish and inter- 
ferometric observations towards 121307 and 122172, aimed 
at investigating their nature. 

In Sect. 13 we describe the observations and the meth- 
ods used for data reduction. In Sect. |31 we present the 
physical parameters of the sources and associated out- 
flows. In Sect.^ we discuss our findings, and give a sum- 
mary and the main conclusions in Sect.|Sl 



2. Observations and data reduction 

The main characteristics of the observed sources are given 
in Table ^ The coordinates corresponding to the sub- 
mm peak position (Molinari et al. 12000(1 in 121307 and 
to the 3.4 mm emission (Molinari et al. I^UU^ in 122172, 
are listed in Cols. 2 and 3; the kinematic distance, d, and 
the Galactocentric distance, i?GCi are given in Cols. 4 and 
5, while the bolometric luminosity, L (from Molinari et al. 
120021 . and the source velocity, wlsr ('f^NHa from Molinari 
et al. I1996|l . are shown in Cols. 5 and 6, respectively. 

Table |2| lists all molecular transitions observed. For 
each line, we give the rest frequency (Col. 2), the tele- 
scopes used (Col. 3), the half power beam width (HPBW) 
(Col. 4), the total bandwidth (Col. 5) and the spectral 
resolution (Col. 6). 



2.1. IRAM-30m observations 

2.1.1. 12C0 (2-1) and C^O (2-1) 

The observations of the ^^CO (2-1) line were carried out 
on February 28 2002 towards 121307, and July 14 2003 
towards 122172. Those of the C^^O (2-1) line were per- 
formed on July 13 2003. We made 5' x 5' maps in the ^^CO 
(2-1) hnc, and 2' x 2' maps in the C^^O (2-1) line using the 
HERA mutibcam receiver in the on-thc-fly mode, in which 
the telescope moves across the source and dumps the data 
with 4" sampling. The dump time was 2 sec per position. 
We used the receiver with a derotator angle of 9.6°. The 
observations were made in position switching mode for the 
^^CO (2-1) line, and in frequency switching mode for the 
C^^O (2-1) line. The reference positions used for the ^^CO 
(2-1) maps were a(J2000) = 2l''55"28!0 and S{32000) = 
48°10'37" for 121307, and a(J2000) = 22''25™55!8 and 
(5(J2000) = 54°35'07" for 122172, and they are CO-free at 
levels of ^ 0.2 K. Telescope HPBW, total bandwidth and 
spectral resolution are listed in Table El Focus and point- 
ing were checked every hour by means of continuum cross 
scans on Uranus or 1253-055. The data were calibrated 
with the "chopper wheel" technique described in Kutner 
& Ulich ( 1981 ). Calibration was checked by measuring the 
^^CO (2-1) line intensity towards the ultracompact Hii re- 
gion W51D. 

2.1.2. CH3C2H, H13C0+ and SiO 

The CH3C2H, H13C0+ and SiO hues were observed to- 
wards 122172 on August 11, 1999. 

For the CH3C2H lines, we simultaneously observed the 
J =(6-5), (8-7) and (13-12) rotational transitions in the 
3 mm, 2 mm and 1.3 mm bands, respectively. The obser- 
vations were made towards the nominal central position of 
122172 given in Tabled We used two spectrometers simul- 
taneously: an autocorrelator, which covered the lowest K 
components with high spectral resolution, and a filterbank 
with low spectral resolution which covered all K compo- 
nents for the (6-5) and (8-7) transitions and up to the 
K = 10 for the (13-12) line. In the following we shall 
make use only of the spectra obtained with the autocor- 
relator because all lines detected fall in the autocorrelator 
band. 

For the Hi3C0+ (1-0) and the SiO (2-1), (3-2) and 
(5-4) {v — 0) lines, 3x3 point maps were obtained, with 
12" (H13C0+) and 24" (SiO) spacing. 

We observed in "wobbler switching" mode, which 
uses a nutating secondary reflector with a beam-throw of 
240" in azimuth and a phase duration of 2 sec. Pointing 
and receiver alignment were checked on nearby point 
sources, and were found to be accurate within 2". The 
data were calibrated with the "chopper wheel" technique. 
Calibration was checked by measuring the intensity of 
Venus at the same frequencies. The main beam brightness 
temperature, Tmb, and the flux density, F^, are related by 
the expression F,j{Jy) = 4.9 Tmb(K). 
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2.1.3. Data reduction and fitting procedure 

CH3C2H is a symmetric-top molecule whose rotational 
levels are described by two quantum numbers: J, associ- 
ated with the total angular momentum, and K, associated 
with its projection on the symmetry axis. This structure 
implies that for each J + 1 — > J rotational transition, 
J + 1 lines are observable (for a detailed description see 
Townes & Schawlow ll975|l . To analyze the CH3C2H data 
we have adopted the fitting procedure described in Fontani 
et al. (|2D02), which assumes the same rest velocity and 
linewidth for all K components. 

All data were reduced with the GAG software devel- 
oped at IRAM and Observatoire de Grenoble. 

2.2. NRAO observations 

Observations of the ^^CO (2-1) transition were made on 
February 21 1998 towards 121307, and on January 10 1998 
towards 122172, using the 12-m telescope of the National 
Radio Astronomy Observatory (NRAO) at Kitt Peak. We 
used the 1 mm SIS array receiver (Payne & Jewell [TMH| . 
consisting of 2 x 4 independent beams. The maps were 
acquired with a 29" spatial sampling, which is approxi- 
mately the telescope HPBW at the ^^CO (2-1) line fre- 
quency. The typical integration time on each position was 
about 4 minutes, with a system temperature of ~ 1000 K 
for 121307, and ~ 2000 K for 122172. Data were taken in 
position switching mode, and the same reference positions 
given in Sect. 12. II were used, which are CO-free down to 
^ 0.2 K. HPBW, spectral resolution and bandwidth are 
hsted in Table El 

2.3. JCMT observations 

The 850/im continuum images were taken on October 16 
1998 with SCUBA at the JCMT (Holland et al. [1998 1 
towards the positions given in Table Q] The standard 64- 
points jiggle map observing mode was used, with a chop 
throw of 2 arcmin in the SE direction. This results in a 
3' X 3' map size. The atmospheric optical depth at 225 
GHz was T225GHZ ~ 0.15. Telescope focus and pointing 
were checked using Uranus and the data were calibrated 
following standard recipes as in the SCUBA User Manual 
(SURF). 

2.4. OVRO Observations 

Millimeter interferometric observations were performed 
using the Owens Valley Radio Observatory millimeter ar- 
ray. 121307 and 122172 were observed in the period from 
October to December 2001. The six 10.4-m dishes were 
employed in three separate configurations (C, L, and H), 
offering baselines from the shadowing limit to approxi- 
mately 220 m. The flexible digital correlator was con- 
figured to observe simultaneously the CO(l-O) and the 
C'^^O(l-O) lines at 2.6 mm. The resolution used was 
^^0.9 km s~^ over a 80 km s^^ bandwidth for the C0(1- 



0), and ^0.65 km s~^ over a 20 km s~^ bandwidth for the 
C^^O(1~0). Continuum observations have been performed 
making use of the 2 GHz bandwidth of the analog correla- 
tor. The observing cycle alternated observations of the two 
target sources and of the phase calibrators (2201+508 and 
BL Lac); a calibration scan was acquired every ^12 min. 
Pointing was checked and refined by frequent pointing ses- 
sions on BL Lac. Bandpass calibration was performed us- 
ing observations of 3C273 and/or 3C345, while the flux 
density scale was determined using observations of cither 
Uranus or Neptune (the expected accuracy is within 15%). 
The raw visibility data were calibrated using the OVRO 
in- house MMA package (Scoville et al. '1993); the cali- 
brated source data were then exported in FITS format 
for further analysis with standard packages. 

2.5. Near-Infrared observations 

2.5.1. Palomar 60-inch broad-band observations 

121307 and 122172 were observed in the three standard 
J, H, and Kg near-infrared broad-band filters using the 
P60IRC (Murphy et al. 1995) camera at the Palomar 60- 
inch telescope on 1999, November 30. For each source we 
obtained a set of dithered frames. After fiat fielding and 
sky subtraction the frames were registered and coadded, 
producing final mosaics of a region of ~3.'5x3.'5 around 
each of the two IRAS sources. Photometric calibration 
was ensured by observations of a set of stars from the 
Arnica and Las Campanas lists (Hunt et al. 119981 Persson 
et al. 11998(1 . The limiting magnitudes of our observations 
were found to be 18.0, 17.5 and 17.0 in J, H, and Kg respec- 
tively. Astrometric calibration was performed using stars 
from the HST Guide Star Catalogue, and later checked 
with 2MASS^ objects in the field (the derived accuracy is 
within 1" rms). 

2.5.2. TNG/NICS narrow-band observations 

Near-Infrared narrow-band observations of 121307 and 
122172 were performed on August 21, 2002, at the 3.56- 
m Telescopio Nazionale Galileo (TNG) at the Roque de 
Los Muchachos Observatory on the Spanish island of La 
Palma. The instrument used is the Near Infrared Camera 
and Spectrograph (NIGS, Baffa et al. l2001|l . a cryogenic fo- 
cal reducer designed and built by the Arcetri Observatory 
IR group as a common-user instrument for the TNG. 
NICS is equipped with a RockweU 1024^ HAWAII near 
infrared array detector. The plate scale used for our ob- 
servations is 0."252/pix, for a total field of view of ^A. 
'5. Series of dithered frames were obtained in the Fell, 
He, H2, and Kc filters, centered at the rest frequencies 
of the [Fell] A ~1.64 ^m transition, the H2(1-0)S(1) 
A ~2.12 /im transition and the nearby continua (see 
Ghinassi et al. 120021 for more details on the NICS filters). 
The seeing at the time of the observations was ~0."9, but 



http://www.ipac.caltech.edu/2mass/ 
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the sky conditions were not photometric, which did not 
aUow to obtain a flux cahbration for our images. 

Data reduction and analysis were performed using the 
IRAF software package, foUowing standard flat-fielding 
and sky subtraction, the individual images were regis- 
tered and the flnal mosaic was produced. The narrow-band 
"continuum" filters were used to subtract continuum emis- 
sion from the Fell and H2 images, assuming that a set of 
stars within each frame were not affected by significant 
line emission or absorption. Note that, because the He fil- 
ter has a central wavelength slightly shorter (^1.57 /im) 
than that of the [Fell] filter, while Kc has a central wave- 
length slightly longer (^2.19 /im) than the H2 filter, very 
red stars will leave a positive residual in the [Fell] contin- 
uum subtracted image and a negative residual in the H2 
continuum subtracted image. Accurate (< O'.'S) astrome- 
try was derived for all mosaics using stellar positions from 
the 2MASS observations. 

3. Results 

3.1. Line emission 

In Figs. H] and we show the ^^CO (2-1), C^^O (2- 
1) and H13C0+ (1-0) (not observed towards 121307) 
spectra taken at the central position in the maps. We 
also show the CH3C2H (6-5), (8-7) and (13-12) spec- 
tra obtained towards 122172 (Fig.OJ: for these lines, only 
single-pointing spectra are available. The ^^CO (2-1) lines 
show non-Gaussian wings, while the C^^O (2-1) line, the 
H13C0+ (1-0) fine and the CH3C2H fines are weU fitted 
by Gaussians. The SiO (2-1), (3-2) and (5-4) {v = 0) 
lines are undetected up to a 3cr level of 0.18 K, 0.4 K and 
0.3 K, respectively. 

Finally, Fig. 01 presents the spectra of the ^^CO (1-0) 
line observed with the OVRO interferometer: both spec- 
tra have been obtained by integrating the emission of the 
line above the 3a level. In this case, the brightness tem- 
perature in the synthesized beam, Tsb, is plotted against 
the velocity. Both spectra clearly show intense emission in 
the wings, and the central peak is lacking: this is due to 
the fact that the channels near the line peak are repre- 
sentative of the more extended emission, which has been 
resolved out by the interferometer. 

The integrated maps of the ^^CO (2-1), C^^O (2-1) 
and H^'^CO"'" (1-0) lines are presented in Figs. Eland |H1 
In 121307, the C^^O integrated intensity map peaks right 
on the sub-mm source detected by Molinari et al. ©OO). 
In 122172 all fine maps are offset by ~ 10" E from the 
millimeter core (Mohnari et al. 2002^. Since the ll^^CO+ 
and C^^O observations have been carried out in different 
observing runs, it is unlikely that this effect is due to a 
pointing error. Moreover, we will see in Sect. IX5TTI that a 
similar offset between line and continuum is confirmed by 
the interferometric observations. 

In Table|21we give the outfiow parameters, namely: the 
velocity range of the blue- and red- wings, Avi, and Av^ 
(Cols. 3 and 4, respectively), the integrated intensity in 
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^ 
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Fig. 1. From top to bottom: spectra of the ^^CO (2-1) 
line towards 121307, observed with the NRA0-12m and 
the IRAM-30m telescopes, and of the C^^O (2-1) fine ob- 
served with the IRAM-30m telescope. All have been taken 
at the peak position in the maps. 

the wings, / Tdv (blue) and (red) (Cols. 5 and 6), and the 
maximum velocity of the wings, Umaxb a-nd Wmaxr (Cols. 7 
and 8), defined as the difference between the maximum 
velocity of the blue- and red- wings and the line peak. 
The latter has been derived from the rest velocity given 
in Table ^ The velocity ranges of the wings have been 
determined from the spectrum of the ^^CO (1-0) line ob- 
served with OVRO, which has the best angular resolution, 
and they are consistent with the wings observed in the 
other lines. The values listed in Cols. 5 and 6 are mean 
brightness temperatures over the blue and the red lobes, 
respectively. 

3.2. Continuum maps 

In Fig. [3 the 850 fim map observed with SCUBA towards 
121307 is superimposed on the 3 mm and the 7 /Ltm emis- 
sion observed with OVRO and ISOCAM, respectively. The 
angular resolution of both the ISOCAM and OVRO maps 
is ~ 6", identical to that of the ^^CO (1-0) observations 
(Table El . The intensity profile of the SCUBA map indi- 
cates the presence of an unresolved central core, and a halo 
extending over ~ 40". Such a halo is not detected in the 
line maps presented in Sect. 13. ll However, it is well known 
that the optically thin emission of dust and the molecular 
line intensities do not depend on temperature and column 
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Table 1. Observed sources: the distance, d, and the bolometric luminosity, L, are taken from Mohnari et ah P2( JU2fl . 
The source velocity, wlsr, corresponds to the molecular gas velocity, derived from NH3 emission (Molinari et al. il996|l . 



Name 


R.A.(J2000) 


Dec.(J2000) 


d 
(kpc) 


Rgc 
(kpc) 


L 

(xIO^Lq) 


«LSR 

(kms-1) 


121307 
122127 


2l''32'"30.6^ 
22*^ 19^08.6^ 


51°02'16.5" 
56°05'02.0" 


3.2 

2.4 


8.4 

7.8 


4.0 
1.8 


-46.7 
-43.8 



Table 2. Observed molecular transitions. 



Molecular 


Rest freq. 


Telescope 


HPBW 


Bandwidth 


Resolution 


transition 


(GHz) 




(") 


(MHz) 


(MHz) 


"^CO (1-0) 


115.271 


OVRO 


6.95x6.12 


26.3 


0.33 


I'CO (2-1) 


230.538 


30m 


11 


320 


0.32 


i^CO (2-1) 


230.538 


NRAO 


29 


150 


0.78 


C^^O (2-1) 


219.560 


30m 


12 


140 


0.078 


Hi3C0+ (1-0) 


86.754 


30m 


29 


70 


0.078 


CH3C2H (6-5) 


102.547(*) 


30m 


24 


70 


0.078 


CH3C2H (8-7) 


136.728(*) 


30m 


18 


70 


0.078 


CH3C2H (13-12) 


222.166(*) 


30m 


11 


140 


0.078 


SiO (v = 0,J=2-1) 


86.846 


30m 


29 


70 


0.078 


SiO {v = 0,J=3-2) 


130.269 


30m 


20 


70 


0.078 


SiO [v = 0,J=5-4) 


217.105 


30m 


12 


140 


0.078 



(*) rest frequency of the K = Q transition 

Table 3. Line parameters for outflow calculations. 



Source 


Line 


Avb 


A«r 


J TAv (blue) 


/■ 


TAv (red) 


^maxj-, 


^maxr 






(km s"^ ) 


(km s"^ ) 


(K 


km s 


) 


(K 


km s~ 


M 


(km s-^ ) 


(km s'^ ) 


121307 


^"CO (2-1) (NRAO) 


(-55, -50) 


(-44, -38) 




10.2 






8.9 




8.3 


8.7 




i^CO (2-1) (30m) 


(-55, -50) 


(-44, -38) 




16.2 






14.9 




8.3 


8.7 




i^CO (1-0) (OVRO) 


(-55, -50) 


(-44, -38) 




13.2 






34.8 




8.3 


8.7 


122172 


I'CO (2-1) (NRAO) 


(-51, -45) 


(-38, -30) 




21.1 






22.5 




7.4 


13.6 




i^CO (2-1) (30m) 


(-51, -45) 


(-38, -30) 




25.4 






30.0 




7.4 


13.6 




'"CO (1-0) (OVRO) 


(-51, -45) 


(-38, -30) 




37.2 






38.1 




7.4 


13.6 



density in the same manner. Furthermore, chemistry and 
optical depth effects can affect the line emission and cause 
a different shape of the maps. 

The OVRO map indicates the presence of a compact 
core, whose diameter is ~ 0.09 pc, located near the emis- 
sion peak of the 850 /im map. One sees that the emission 
peak of the fSOCAM map is offset by -- 6" in R.A. from 
the the 3 mnr and the 850 /xm emission peaks. This is con- 
firmed by the 15 /im ISOCAM map, whose emission peak 
is nearly coincident with that detected at 7 /im. Therefore, 
the ISOCAM maps may suggest that in 121307 the bulk of 
the MIR emission does not arise from the millimeter core. 
However, we point out that these maps have to be re- 
garded with caution because the astrometric uncertainty 
may be as large as ~ 10" (see Cesarsky & Blommaert 
I2000|l . and maps with better astrometric precision and 
angular resolution are required to assess the origin of the 
mid-infrared emission. 



with the MIR emission, we have superimposed the MSX^ 
21 /im image, because no ISOCAM maps are available for 
this source. The bottom panel presents an enlargement of 
the central region with a contour plot of the 3 mm con- 
tinuum observed with OVRO, which shows the presence 
of a compact millimeter core (~ 0.04 pc) approximately 
located at the center of the SCUBA map. Although the 
angular resolution of the MSX map is '^ 4 times worse 
than that of ISOCAM, we find that also in this case the 
peak of the MIR emission does not coincide with that at 
850 /im and 3 mm. 

Near-infrared images of both sources are presented in 
Fig. 1^ and 1101 obtained at 2.2 /im (Kc-band) and 1.7 /xm 
(He-band). In both maps, the half maximum power con- 
tour of the 3 mm core, and the MIR emission are also 
shown. The Kc-band image of 121307 (top panel of Fig.O 
shows a cluster of stars, the "centroid" of which is located 
at ~ 6" (~ 0.08 pc) in R.A. from the peak of the millime- 



The top panel of Fig. |H| shows the SCUBA map ob- ^ ^g^ images have been taken from the on-line MS X 
served towards 122172. In this case, to make a comparison database http://www.ipac.caltech.edu/ipac/msx/m8x.htmll 
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Table 4. Physical parameters of the 121307 core. 9s and D are the angular and linear diameter, respectively, computed 
for a distance of 3.2 kpc. 





Os 


D 


Mvir 


Mh, 


nH2 




(arcsec) 
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Table 5. Physical parameters of the 122172 core. Og and D are computed for a distance of 2.4 kpc. 
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ter core. However, the brightest cluster member in Kc is 
~ 1" South of the OVRO peak. The 15 /im ISOCAM peak 
lies on the western side of the cluster. It seems probable 
to us that such emission is due to a combination of emis- 
sion from the cluster members and small particle emission 
from the adjacent cloud. 

In 122172, the 3 mm core is nearly at the center of 
a cluster of NIR sources distributed as a patchy "ring" , 
detected both at 1.7 and 2.2 /im (Fig. I10() . In this case 
also the 21 /im emission is roughly coincident with the 
position of the cluster. 

It is interesting to notice that in both sources the stars 
detected in the near-infrared images are located close to 
the millimeter core, but have a distribution that does not 
exactly coincide with it: the "average position" of the red- 
dest stars of the cluster is ~ 6 — 8" offset from the mil- 
limeter peak, i.e. well above the positional uncertainties 
on the maps, estimated to be < 1". This is reminescent of 
the observations of the candidate massive protostar IRAS 
23385+6053 (Fontani et al. ^^, in which a compact 
massive core, believed to host a newly formed massive 
(proto)star, is surrounded by a cluster of more evolved B 
stars, distributed as a patchy ring around the core. 

3.2.1. The NIR clusters 

The near infrared broad band images of 121307 and 122172 
clearly show a concentration of red star-like sources near 
the IRAS positions. This is consistent with the detection 
of young stellar objects closely associated with the molec- 
ular clumps. Following the methods described in Testi et 
al. IJ19(T71 I1998|l we have verified the presence of a stel- 
lar surface density peak and we have estimated the rich- 



ness indicator Iq for both regions. The resulting values 
are Ic ~13±3 for 121307 and Ic -^9±3 for 122172. 

In order to compare these numbers with other cluster 
surveys, it is necessary to give a rough estimate of the 
minimum stellar mass we are able to detect in our obser- 
vations. Following Testi et al. H1998|l . but using the Palla 
& Stabler H1999|l pre-main sequence evolutionary tracks, 
we have computed the expected stellar mass sensitivity of 
our observations. Neglecting the effects of infrared excess 
and assuming an age of ~ 0.1-1 Myr and an extinction 
Av'=10 mag, we estimate that our observations should 
be sensitive to young stellar objects more massive than 
~0.5-1.5 Mq. This is necessarily a rough estimate, and 
our observations are certainly not complete to these low 
masses. 



3.3. SED and luminosities 

In Figs. ^] and E| we plot the Spectral Energy 
Distribution (SED) of both sources. 

It is clear that these regions are very complex and a 
proper modeling of the SED would require a detailed ra- 
diation transfer calculation, taking into account both the 
geometry of the core and cluster as well as the radiation 
field of each cluster member and the contribution of very 
small dust grains and PAHs. This approach would require 
a substantial number of assumptions and goes beyond the 
scope of this paper. Our aim here is to estimate the rel- 
ative contribution to the SED of the two main compo- 
nents detected in each region: the NIR cluster and the 
cold molecular core. This simplified approach is described 
in the following. 
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Fig. 2. Spectra of the -^^CO (2-1) line observed with the 
IRAM-30m telescope and the NRAO telescope, and of the 
H13C0+ (1-0) and C^^O (2-1) lines observed with the 
IRAM-SOm telescope. All spectra refer to the peak posi- 
tion in the maps. 



Let us consider first 121307: the SED of this source 
was already derived by Molinari et al. H2U0U|I . With re- 
spect to their Fig. 3, Fig. [TTl contains more data, namely 
MSX, IRAS, OVRO and VLA. On the basis of the ob- 
servations presented in Sect. 13.21 we can distinguish be- 
tween two components: the 3 mm core (hereafter 121307- 
C), and the emission of a stellar cluster, detected in the 
NIR images. The millimeter and sub-millimeter flux den- 
sities arise from the core I21307-C, while, as discussed in 
Sect. 13.21 the MIR data are likely due to the nearby clus- 
ter. Therefore in Fig. llll we have performed two grey-body 
fits. 

The solid line represents the best grey-body fit ob- 
tained using data representative of the emission arising 
from I21307-C, i.e. from millimeter to far-infrared wave- 
lengths. Assuming a dust opacity index f3 — 2, which is 
a value typically found in envelopes of high-mass YSOs 
(Hunter 1997; Molinari et al.[2D003 the corresponding pa- 
rameters are: dust temperature T = 30 K, angular diam- 
eter = 10", and mass M ~ 90 Mq. This fit yields a 




V|SR (km s 



Fig. 3. From top to bottom: spectra of the CH3C2H (6- 
5), (8-7) and (13-12) lines observed towards 121307 with 
the IRAM-30m telescope. The numbers under the spectra 
indicate the position of the different K components. The 
Vlsr is relative to the K = line. 



luminosity of 4 x IO'^Lq, which corresponds to the lumi- 
nosity of a B1.5 ZAMS star. This luminosity has to be 
regarded as an upper limit to the "real" luminosity of the 
core I21307-C. This is because the IRAS flux densities at 
100 and 60 yum may arise both from the core and from the 
external cluster, because the IRAS beam is much larger 
than the core region. However, the goodness of the grey- 
body fit shown in Fig.llllsuggests that the 100 and 60 /um 
flux densities are mainly due to I21307-C, and therefore 
the luminosity computed above is close to the "real" core 
luminosity. 

The dotted line is a grey-body fit to the data at wave- 
lengths A < 25 iim, and it represents the emission coming 
from the external cluster. The fit is obtained for a temper- 
ature of 150 K. By integrating the flux densities we obtain 
a luminosity of 3.6 x lO'^L©. 

Let us now consider 122172: Fig. IT^ shows the SED of 
this source. The images obtained at different wavelengths 
show that also in this case there is a molecular core (here- 
after I22172-C), detected in the millimeter continuum and 
molecular lines maps, and a cluster of stars revealed in the 
near-infrared images. Therefore, as for I21307-C, we have 
performed a two grey-body fit to the SED. Using f3 — 2, 
from Fig. El we have derived the following best fit param- 
eters: dust temperature T = 27 K, angular diameter of the 
source = 9" and mass of ^ 83 M©. The resulting luminos- 
ity is 2.2 X IO'^Lq. In this case also, this luminosity must be 
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Fig. 4. Top panel: spectrum of the ^^CO (1-0) line ob- 
served with the OVRO interferometer towards 121307. 
The spectrum has been obtained integrating the emission 
over the 3(t level. Bottom panel: same as top panel for 
122172. 



regarded as an upper limit for the luminosity of I22172-C, 
but, as discussed for 121307, the fact that the fit is reason- 
able suggests that the emission at 100 and 60 /im arises 
mainly from the core, so that the real core luminosity is 
close to the value estimated before. The dotted line shows 
the grey-body fit to data with A < 25 /im, which represent 
the emission due to the cluster of stars surrounding the 
core. It has been obtained for a temperature of 140 K and 
a luminosity of 2.2 x lO^i©. 

The shape of both spectra is similar to that of the 
massive protostar candidate IRAS 23385+6053 (Fontani 
et al. EI 



1^ . In addition, as pointed out in Sect. 13.21 the 
source structures are also similar. This may suggest that 
the evolutionary stage of the sources is the same. We will 
come back to this point in Sect. 14.21 

3.4. Mass, density, and temperature estimates 

Core masses have been derived in two different ways: from 
the continuum emission and from the lines. In the first 
case, assuming dust emission to be optically thin, we can 
compute the core mass M directly from the observed flux 
density: 



M = 



F,d^ 
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Fig. 5. a: Map of the ^^CO (2-1) Hne integrated between 
—55 and —38 km s~^, observed towards 121307 obtained 
with the NRA0-12m telescope. Contour levels range from 
60 to 200 by 30 K km s~^. The thick contour indicates 
the half maximum power contour. The cross at the map 
center shows the position of the sub-mm peak (Molinari et 
al. l2000|l . b: Same as a for the IRAM-30m telescope. Same 
contour levels are shown, c: Same as a for the C^^O (2-1) 
line integrated between —50 and —42 kms~^, observed 
with the IRAM-30m telescope. Levels range from 2.2 to 
9.2 by 1.0 km s'^ 



where F^ is the flux density integrated over the area de- 
scribed by the 3cr level of the SCUBA map, d is the 
distance, k^ is the dust opacity parametrised as k^ = 
/c^o (I//230 GHz)'^ (in units of cm^ g~i), B{v,T) is the 
Planck function and T is the dust temperature. We as- 
sume fcjy — 0.005 cm^ g^^ (Kramer et al. I1998|l . that im- 
plies a standard gas-to-dust ratio of 100. We have also 
deduced the core mass by means of the total column den- 
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Fig. 6. a: Map of the ^-^CO (2-1) line integrated between 
—51 and —30 km s"-"^, observed towards 122172 with the 
NRA0-12m telescope. Contour levels range from 80 to 
255 by 35 K km s~^. The thick contour indicates the half 
maximum power level. The cross at the map center shows 
the position of the 3.4 mm peak (Molinari et al. l2002|l . b: 
Same as a for the IRAM-30m telescope. The same con- 
tour levels are shown, c: Same as a for the H^^CO"'"(1-0) 
line, observed with the IRAM-30m telescope, integrated 
between —46 and —40 km s^^. Contour levels range from 
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Fig. 7. SCUBA map of 121307 at 850 /im (grey-scale) su- 
perimposed to the 3 mm OVRO map (solid lines) and the 
7 ^m ISOCAM map (dashed lines). Contour levels range 
from: 0.2 (4cr) to 1.1 by 0.16 Jy beam"! in the SCUBA 
map; from 0.003 (- 3cr) to 0.015 by 0.001 Jy beam^i in 
the OVRO map; from 20 to 212 by 24 Jy beam^i in the 
ISOCAM map. 



sity derived from some of the lines observed, according to 
the expression: 



Mine — 



4 X 



(2) 



where D is the source diameter, mn^ the mass of the 
H2 molecule, X the abundance of the molecule relative 
to H2, and A'tot is the source averaged column density. 
For the C^^O molecule, following Wilson & Rood l|1994|l 
we assumed that the oxygen isotope ratio ^^O/^^O de- 
pends on the galactocentric distance Rqc^ according to 
the relationship i^Q/iSO = 58.8i?Gc(kpc) + 37.1; using 
the galactocentric distances given in Tabled and assum- 
ing a mean abundance of 10^^ for the ^^CO molecule, 
we infer a C^*0 mean abundance of ~ 1.9 x 10~^. For the 
H^'^CO^ molecule, we applied a similar argument: the car- 
bon isotope ratio, ^'^C/^^C, depends on the galactocentric 
distance Rgc, according to the relationship ^^C/^'^C — 
7.5i?Gc(kpc) + 7.6 ; thus, assuming a mean abundance of 
X = 10-9 for the HCO+ molecule (Irvine et al. lTMTjI . we 
derive a H^^CO+ mean abundance of ^ 1.5 x 10"^^. 

Furthermore, from the linewidth and the core linear di- 
ameter of each molecular line, we can estimate the virial 
mass: assuming the source to be spherical and homoge- 
neous, neglecting contributions from nragnctic field and 
surface pressure, the virial mass is given by: 



AfviR(Af0) 



0.509 d(kpc) Os (arcsec) Aui/2(kni/s) 



(3) 



The mass estimates deduced from Eqs. ^ and ^ are 
listed in Col. 5 of Tables 0] and [S] while the virial masses 
are shown in Col. 4. In Col. 6 we also give the mean H2 vol- 
ume density, nH2, deduced from the corresponding mass 
and diameter. The values found clearly show that both 
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1.5 by 0.24 Jy beam~^. Bottom panel: enlargement of the 
SCUBA map around the central region. Contours indicate 
the OVRO 3 mm emission, whose levels range from 0.004 
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Fig. 9. Kc-band (top panel) and He-band (bottom panel) 
images of 121307. The position of the 3 mm core is indi- 
cated by the solid line, which represents the half maxi- 
mum power contour of the OVRO maps shown in Fig. {7\ 
The dashed lines indicate the ISOCAM 15 /im image. The 
same contour levels as in Fig. [3 are shown. 



cores are massive, and the masses needed for virial equilib- 
rium are within a factor of ~ 2 from those computed from 
the continuum emission or from the lines. Hence the cores 
are likely virialised. Furthermore, nu^ increases towards 
the core center in both sources; the increment is roughly 
consistent with what is expected for a density profile of 
the type riHs oc R~^, which is the typical profile found in 
molecular clumps that surround high-mass YSOs (see e.g. 
Fontani et al. 2002, Beuther et al. l2002b|l . As usual, the 
uncertainties on the mass estimates are mainly due to the 
molecular abundances, the assumed dust opacity and the 
gas-to-dust ratio, which are known to be affected by large 
errors and difficult to quantify. 

One notices that the masses estimated from the flux 
densities at 850 /im arc larger than those derived from the 
C-'^^O (2-1) line, although the source diameters are bigger 
in the C^^O (2-1) line. This may be due to the fact that 
all diameters have been estimated from the FWHM of the 



emission, and in the 850 /im maps a significant fraction of 
the total integrated emission arises from a halo of '^40" 
in diameter, which is not detected in the line maps. 

Finally, from the CH3C2H spectra observed towards 
122172 (Fig. I2I1, we have derived the kinetic temperature 
and the total column density of this source by means of the 
rotational diagram method (see e.g. Fontani et all2 002|l . 
The fundamental assumption of this method is that the 
gas is in local thermodynamic equilibrium (LTE) condi- 
tions, which we believe to be valid for CH3C2H because 
of its low dipole moment. If we further assume the lines 
to be optically thin, one can derive the column density A'i 
of the upper level of each transition from the integrated 
intensity of the line. We obtain the source kinetic temper- 
ature and column density listed in Cols. 1 and 2 of TableEl 
respectively. We also list the virial mass (Col. 3), the gas 
mass (Col. 4) and the H2 volume density (Col. 5), com- 
puted as explained before for the other molecular tracers. 



Fontani et al.: IRAS 21307+5049 and IRAS 22172+5549 



11 



56°05'2G' 



56°05'00' 



O 

o 
o 

CM 



'^ *l ' 


JHh 


__L^ ' 


l-^ ^ 








_ — 


\ 








; ^-'"Z*- 


— -m 




1 .'"Ji— 

1 / r"^ *"-- 


Hf" 


">•:' • - 


/ ' ^~~~ 


\"^ \ 


•/ ' 


'ilfti 






-' '.iWt 






\ \ IT^ 


11/ 






III' 
1 1 \ 


/ 

; * 


,SN^^ *^-'' 


1 1 1 




\ \ V ^ ^X ^ 


^ ; 




.w^^-:-^'' 


y ^1 


/ 12^172 




^^ 


/ Kc 








\, ^4_ , .-' 


, ^i 


, 



22 19 10' 



08 



56°05'20' 



O 
CD 
Q 



56°05'00' 



I [ I \ i i J / / ' ' I / 



\ \ 



He 



22 19 10 



08 



R.A. (.2000) 

Fig. 10. Same as Fig.Elfor 122172. The dashed lines repre- 
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as in Fig. |S1 

Table 6. Temperature, mass and H2 volume density es- 
timates obtained from CH3CCH for 122172. A CH3C2H 
abundance of 2 x 10~^ has been assumed. As source di- 
ameter, we have used that of the SCUBA map. 

^71 1^1 l^r Mh^ n^2 

(K) (cm-2) (Mo) (Mo) (cm-^) 



27 4 10" 



73 



92 



1.0 IQ' 



The diameter used to compute jt-Hs is that derived from 
the SCUBA map. We assumed an abundance of CH3C2H 
relative to H2 of 2 x 10^^. This is an average value taken 
from Fontani ct al. ( 2002) and Wang et al. IJTMSJI . who 
studied this molecule in massive star forming regions. 

3.5. The outflows 

3.5.1. Integrated maps and morphology 

Figures El and E| show the integrated intensity maps of 
the red and the blue wings of the ^^CO (2-1) line. Panels 
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Fig. 11. SED of 121307. The meaning of the symbols is 
given in the Figure. The arrows at 2 mm and 2 cm rep- 
resent the upper limits (3cr level) of the JCMT and VLA 
maps, respectively. The solid line represents a grey-body 
fit to the far-infrared and millimeter points, with dust tem- 
perature r = 30 K and opacity index /3 = 2. The dotted 
line is a grey-body fit to the data at A < 25 /xm, whit dust 
temperature T = 150 K. 
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as in Fig. 1111 The arrows at 2 cm and 6 cm represent 
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opacity index f3 — 2. The dotted line is a grey-body fit to 
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Fig. 13. Molecular outflow observed towards 121307. a: 
i^CO (2-1) line observed with the NRA0-12m telescope, 
integrated under the blue-shifted (solid line) and the 
red-shifted (dashed line) wings. The velocity ranges are 
(-55,-50) kms-i (blue) and (-44,-38) km s^^ (red). 
Contour levels begin from 6 {^Aa) K km s^^, with a spac- 
ing of 2 K km s~^. The overlaid grey-scale represents the 
SCUBA 850 /im emission map: the same contour levels 
as in Fig. [7| are shown, b: Same as a for the ^^CO (2- 
1) line observed with the IRAM-30m telescope. Contour 
levels begin from 2.5 K km s~^ (~ 3cr level in the maps) 
and have a spacing of 2 K km s^^. c: Same as a and for 
the ^^CO (1-0) line observed with OVRO. Contour values 
begin from 0.9 K km s~^ (~ 3cr) and have a spacing of 
0.3 K km s^^. The grey-scale represents the core 121307- 
C observed with OVRO at 3 mm, with the same contours 
as in Fig.[71 



a and b show the observations performed with the NRAO- 
12m telescope and the IRAM-30m telescope, respectively. 
Solid and dashed lines correspond to the integrated emis- 
sion under the blue and the red wings, respectively. The 
velocity intervals are given in Table O These have been 
superimposed on the SCUBA map at 850 /im, which in- 
dicates the position of the cores. 

The outflow associated with 121307 is oriented approx- 
imately NW-SE, with the blue-shifted gas towards the SE, 
and the lobes appear to be well separated. In contrast, to- 
wards 122172 the lobes overlap completely. However, in 
both cases, the angular resolution is too poor to resolve 
the outflow structure. The interferometric observations 
obtained with OVRO in the ^^CO (1-0) hue give us better 
understanding of the structure of the flows: in panels c of 
Fiefs. 1131 and 1141 the outflow maps are shown. These have 
been obtained integrating the line emission over the same 
velocity intervals as in panels a and b. 

Let us consider the two outflows separately. 121307 
(Fig. Et) presents an outflow oriented NW-SE, with the 
blue-shifted gas located to the SE. The lobe peaks are 
centered near the core position. The morphology is quite 
chaotic and the outflow seems to be made at least of two 
components: a stronger one centered on the core 121307- 
C, and a weaker one, offset by ~ 15" to the west. Various 
authors (e.g. Beuther et al. l2002;ajl have demonstrated the 
presence of multiple outflows in massive star forming re- 
gions. Hence, such a complex morphology may be due to 
the fact that we are seeing the superimposition of several 
flows. 

In 122172 (Fig.lMbV the outflow axis is clearly oriented 
N-S, with red-shifted gas to the south and blue-shifted 
gas to the north. The whole outflow has an extension of 
~ 30", corresponding to ~ 0.35 pc at 2.4 kpc. In this 
case the lobes have similar sizes and there is little overlap 
between them. It is striking that the center of the lobes 
is offset by ~ 5" from the dust emission peak. This sug- 
gests that the source driving the outflow is not located at 
the center of the dust emission. The outflow might origi- 
nate from a source belonging to the more evolved cluster, 
and not from the massive object embedded inside 122172- 
C; in fact, comparing the image at 1.7/xm (bottom panel 
of Fig. I10|l with the panel c of Fig. El one notices that 
the center of the lobes roughly corresponds to one of the 
near-infrared sources. Thus, the molecular core might host 
a source in a very early evolutionary stage prior to the for- 
mation of an outflow. However, as we will briefly discuss 
in Sect. 13.5^ the flow parameters indicate that the source 
driving the outflow must be massive. Therefore, wherever 
the source driving the outflow is located, it must be a 
massive young object. Since the stars of the cluster are 
probably low- to intermediate-mass stars, it is unlikely 
that the flow originates in one of the cluster members. 

We can compare our OVRO maps with those presented 
in Molinari et al. l|M7^ . who used the HCO+ (1-0) line 
to observe the outflows in both sources. They adopted 
velocity intervals slightly different from ours. After inte- 
grating the ^^CO (1-0) lines in the same intervals, one 
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sees that our maps are consistent with those of Molinari 
et al. P002) for 121307, but not for 122172. This might 
be explained by the fact that Molinari et al. H2002|l ob- 
served only the inner part of the region observed by us 
using a more compact configuration of the interferometer, 
sensitive to more extended structures. 

3.5.2. Physical parameters 

The characteristics of the flows are listed in Table [7| in 
which the source and the tracer used are given in Col. 1 
and 3, respectively. The derived physical properties have 
been computed using the same relations given in Beuther 
et al l2()()2ai the mean column densities of the blue- and 
the red- lobes, A'b and N,-, are given in Col. 4 and 5, the 
mass of the blue- and the red- lobes Mt, and M^, and the 
total mass, Mout, hi Col. 6, 7 and 8, respectively, the mo- 
mentum Pout in Col. 9, the energy E in Col. 10, the size 
r in Col. 11, the dynamical timescale idyn in Col. 12, the 
mass entrainment rate Mout in Col. 13, the mechanical 
force Pout in Col. 14 and the mechanical luminosity L,n in 
Col. 15. All values have been derived assuming the mean 
inclination angle of 57° (Cabrit & Bertout .1992^) with re- 
spect to the line of sight. 

To determine the size r, we have adopted the maximum 
extension from the center of the lobes (Cabrit & Bertout 
I1992II . Since in 121307 the lobes have different extensions, 
we have adopted their arithmetic mean. 

The velocity intervals chosen for the line wings are the 
same used for the OVRO ^^CO (1-0) spectra (see TableEl). 

The ^^CO (2-1) and (1-0) lines may have large opac- 
ities also in the wings. Hence, we should take into ac- 
count this effect in computing the column density. For the 
^^CO (2-1) line we have followed the method outlined in 
Beuther et al. H20()2all : one can assume a mean value for 
the i^CO/i^CO (2-1) line wing ratio of - 0.1 (Choi et 
al. I1993|l . and that this ratio is constant throughout the 
outflow (Cabrit & Bertout 11986(1 . Even though we have 
not made maps of the ^^CO (2-1) line, we can verify if 
this assumption holds for our sources using the data of 
Brand et al. ^2001) : they observed the ^^CO (2-1) fine 
towards 121307, and by comparing their spectra with ours 
in the ^^CO (2-1) line, we conclude that a mean ratio of 
0.1 is plausible for the wings. Under this assumption, the 
H2 column densities of the wings, N^ and N^, have been 
computed using the relation (see Beuther et al. l2002ia|l : 



N == 



H, 






13C0; 47r3/li^2^2 



exp(-16.6/Te,) 



0.1 



TMB''C0(2-l)dw, 



(4) 



wing 



We have assumed tI^ = 89 x 10* and Tox = 30 K as 
in Beuther et al. H2002a|l . to make a consistent compari- 
son. For this reason the correction for the galactocentric 
distance has not been applied. However, such a correction 
would not affect any of the conclusions drawn by us and 
Beuther et al.. 



For the ^^CO (1-0) line, we have estimated the optical 
depth in the wings as follows: from the mean line wing 
ratio ^^CO/^^CO (2-1)= 0.1, and a relative mean abun- 
dance 



H2 



13 CO 



= 89 X 10 , we derive an optical depth of the 
^^CO (2-1) line in the wings of r ~ 5. One can demon- 
strate that, in LTE conditions, for a kinetic temperature 
of 30 K, the optical depth of the (2-1) line is ~3 times 
that of the (1-0) line. Hence, we adopted a mean value 
T = 1.7 for the ^^CO (1-0) line to compute the column 
density in the wings, and corrected the column density 
estimates accordingly. 

The parameters deduced from the ^^CO (1-0) line 
are in good agreement with those previously found by 
Molinari et al. (jMSJ, who observed the HCO+ (1-0) line 
towards both sources with similar angular resolution, as 
mentioned in Sect. 13. 5m The outflow masses, the mechan- 
ical forces and the mechanical luminosities are at least an 
order of magnitude greater than those typically associ- 
ated with low-mass YSOs (see e.g. Bontemps et al. 11996(1 . 
and in good agreement with those found by Beuther et 
al. (2002aJ in a sample of massive YSOs, believed to be 
precursors of UC Hll regions. A detailed comparison be- 
tween our results and those found by other authors (both 
in low-mass and high-mass star forming regions) will be 
done in Sect. 14. ll 



3.5.3. Near- Infra red H2 and [Fell] images. 

Fig- El presents the image of 122172 in the H2 and [FeH] 
lines. They have been superimposed on the outflow image 
derived from the ^^CO (1-0) line. Also indicated in Fig. 1151 
is the position of the core I22172-C. The [FeH] emission 
is in the center of the blue lobe of the outflow, and it 
is oriented approximately NE-SW. The H2 emission has 
roughly the same inclination, and it seems to follow the 
boundary of the [FeH] emission. Hence, we conclude that 
both these tracers are associated with the ^^CO flow seen 
with OVRO. 

It is well known that NIR lines of H2 and [Fell] are 
tracers of the hotter and faster component of a flow in low- 
mass YSOs (e.g. Nisini et al. 12001(1 . In particular, [Fell] is 
thought to be a tracer of the inner region of the jet, while 
the H2 emission is likely due to the shocked gas on the 
surface between the jet and the ambient material. From 
Fig. El one can see that our observations are consistent 
with this picture: the H2 and [Fell] emission both arises 
from the inner parts of the outflow lobes. 

No H2 or [Fell] line emission was detected in 121307. 



4. Discussion 

In order to understand the nature of IRAS 21307+5049 
and IRAS 22172+5549, in the following we wiU discuss 
the results obtained in Sect. 13 
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Table 7. Physical parameters of the outflows. Column densities of the blue- and the red- wing (A'b and N,-) are 
expressed in units of lO^^cm"^, masses (Mb, M^ and Mout) in Mq, momentum (Pout) in M0 km s~^, energy (E) in 
10**^ ergs, size (r) in pc, dynamical timescale (^dyn) in 10^ yr, mass entrainment rate (Mout) in IO^^Mq yr~^, mechanical 
force (Pout) in lO~'^M0km s~^ yr~^ and mechanical hmiinosity (Lm) in Lq. All values have been calculated assuming 
the mean inclination angle of 57°. 



Source 


tracer 


iVb 


K 


A4 


M, 


Mout 


"out 


E 


r 


^dyn 


Mout 


-Pout 


J-^m 


121307 


'^CO (1-0) (OVRO) 


33.8 


17.9 


2.3 


0.8 


3.2 


48.3 


3.7 


0.1 


1.3 


4.5 


2.9 


10.5 




12 CO (2-1) (30m) 


0.8 


0.7 


5.2 


3.0 


8.2 


127.8 


20.0 


0.7 


4.3 


1.9 


2.9 


6.8 




i^CO (2-l)(12m) 


0.5 


0.4 


15.8 


4.2 


20.0 


311.0 


48.3 


1.0 


5.7 


3.4 


5.5 


13.1 


122172 


I'CO (1-0) (OVRO) 


52.4 


39.1 


1.1 


1.6 


2.7 


55.6 


12.4 


0.1 


0.7 


7.6 


12.2 


14.2 




12 CO (2-1) (30m) 


1.2 


1.6 


5.0 


3.4 


8.4 


153.5 


31.0 


0.4 


1.7 


4.7 


9.0 


27.3 




'^CO (2-l)(12m) 


1.0 


1.1 


12.8 


10.1 


22.9 


428.7 


89.3 


0.6 


2.9 


7.8 


14.8 


79.3 



4.1. Comparison with outflows associated with other 
YSOs 

In Fig. ^1 we plot the mechanical force, Pouti of the 
outflows against the bolometric luminosities Pboi of the 
sources. Our data are compared to those derived by 
Bontemps et al. II1996|I . who studied outflows associated 
with low- mass objects, and those of Beuther et al. I)2()()2a|l . 
who analysed massive molecular outflows. We also com- 
pare our results to those obtained by Molinari et al. 
H1998b|l and Fontani et al. iPTOjl for IRAS 23385+6053, a 
candidate massive protostar. Figure lTBl is similar to Fig. 4b 
of Beuther et al. (f2002al : however, unlike these authors, 
we consider only sources without ambiguity between the 
near and the far kinematical distance. Also the values com- 
puted by Beuther et al. (2002a) have been corrected for 
the mean inclination angle of 57°: this allow us to make a 
consistent comparison. 

Beuther et al. I)2002a|l have shown that there is a cor- 
relation between the mechanical force Pout and the lu- 
minosity of the core driving the outflow. They suggest 
that a similar relation, established for low-mass objects 
(Bontemps et al. ll99(jf) . holds up to the high-mass regime. 
Our results and those obtained by Fontani et al. H2004|l 
for IRAS 23385+6053 appear to confirm this statement. 
In particular, 121307 and 122172 have bolometric lumi- 
nosities close to those of early B stars. 

Finally, we have computed the dynamical timescale of 
the outflows. These values represent a lower limit to the 
outflow age (see e.g. Parker [TMT|l . For 121307 and 122172, 
we have found (see Table E)) that idyn is of the order of 
10'' — 10^ yr. Beuther et al. H2()()2a|l found values ranging 
from 2 X lO'* to 2 x 10^ yr. These values are again in good 
agreement with our estimates. 

4.1.1. Mass accretion rates and accretion efficiencies 

From the mass entrainment rate, we can derive the mass 
accretion rate, Maco and the accretion efflciency, /ace- 
For both, we will adopt the same argument outlined by 
Beuther et al. (|2002a|l . 

Based on the assumptions that the source driving the 
outflow is a single star, the momentum of the outflow is 



equal to that of the internal jet and that the ratio be- 
tween mass loss rate of the jet and mass accretion rate 
is Mjet ^ 0.3Afacc (Tomisaka UnnHl Shu et al. HM^ . we 
obtain accretion rates of order lO~^M0yr~-' (see Table|Hl). 
Such values are slightly higher than those derived by 
Beuther et al. H2002a|l in sources with comparable lumi- 
nosities (from ~ 10^ to ^ 10^ Lq), but still consistent 
within the uncertainties. 

From the accretion rate we estimate the accretion ef- 
ficiency, /ace: 



/acc — 



Ma 



Mcoro/te 



(5) 



Using Mace, iff, and the values of Afcorc computed from 
the 850 ^m map, we derive /acc ~ 0.03 and ~ 0.08 for 
121307 and 122172, respectively (see Table 0. 

/acc can also be determined following the argument 
outlined in Richer et al. 1)2000(1 . They consider two ap- 
proaches: 

(a) the luminosity is accretion dominated: 

L Wjot 



/acc = 2 X 10- 



Poutfkep ^'kep 

(b) the luminosity is mainly due to ZAMS stars: I2002a|l 

— 3 -''^*''^kcp 



(6) 



u 



2 X 10" 



(7) 



where Wkcp — \/GM^/R^ is the stellar escape velocity, 
assumed equal to the jet velocity, and M* and P* are the 
stellar mass and radius for a ZAMS star with bolometric 
luminosity L. The luminosities used are those derived in 
Sect. 13.31 Pout and Pout have been calculated as explained 
in Sect. 1^3^ 

The values of /acc computed from Eqs. (jSJ and lO 
are listed in Table |H1 together with the result for IRAS 
23385+6053 and the average values deduced for the 
sources of the Beuther et al. sample with comparable lu- 
minosity. 

4.2. Nature of the sources 

We stress that the quantities previously obtained have 
been derived under several hypotheses, the most impor- 
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tracer 



Mace 

(10-* Meyr-i) 



/ace ( ) /ace (a) /aec (b) 



121307 




i^CO (1-0) 


0.8 


0.04 


0.01 


0.4 






i^CO (2-1) (30-m) 


0.6 


0.03 


0.01 


0.2 






i^CO (2-1) (12-m) 


0.4 


0.02 


0.02 


0.2 


122172 




i^CO (1-0) 


1.3 


0.07 


0.003 


0.3 






i^CO (2-1) (30-m) 


1.4 


0.08 


0.003 


0.1 






I'CO (2-1) (12-m) 


0.8 


0.05 


0.006 


0.09 


IRAS 23385+6053 (*) 


HCO+ (1-0) 


1.5 


0.006 


0.01 


0.1 


Beuther et al. 


(*) 


i^CO (2-1) (30-m) 


0.3 


0.02 


0.08 


0.2 



(*) Molinari et al. Ill998bll 

( ) average values for the sources of the Beuther et al. 12002at sample with luminosity from 10'^ to 10* L© 

(*) computed from Eq. (|SJ| 

(a) computed from Eq. © 

(b) computed from Eq. (|7J 



tant being the assumed jet velocity and the efficiency of 
the ejection mechanism (i.e. the ratio between mass loss 
rate and accretion rate, assumed to be 0.3), which are both 
very uncertain. Hence, the values listed in Table|Hlhave to 
be taken with caution. Neverthless, the results presented 
in Sects. EU and UTTl indicate that 121307 and 122172 
have physical parameters comparable to those obtained 
by Beuther et al. H2002a|l in sources with similar lumi- 
nosities, and hence likely in similar evolutionary stages. 

In Fig. El we plot /ace estimated using Eq. ((SJ , ® 
and for 121307, 122172 and IRAS 23385+6053 (big 
symbols), and for the sources of the Beuther et al. sam- 
ple with comparable luminosity (small symbols) . Although 
the spread of the values is about an order of magnitude for 
the three estimates of /ace, one can notice that the values 
derived from Eq. Q are more consistent with those com- 
puted from Eq. ©, i.e. under the assumption of accretion 
dominated luminosity, than with those estimated in the 
case of luminosity dominated by a ZAMS star. This result 
favours the conclusion that our sources, as well as those 
of Beuther et al. , are likely massive protostars still in the 
accretion phase. 



4.2.1. Comparison with the massive protostar 
candidate IRAS 23385+6053 

In Sects. 13.21 and 13.31 we have shown that the continuum 
maps indicate the presence of two "sub-regions" in both 
121307 and 122172: a massive compact core (called 121307- 
C and I22172-C) and a cluster of stars which surround the 
core. From the maps and the SEDs, we concluded that the 
core is responsible for the continuum emission at wave- 
lengths longer than ^ 25 /im, while the emission at shorter 
wavelengths is mainly due to the cluster. Grey-body fits 
to the SEDs indicate that the cores have luminosities of 
^ 10^ Lq, an order of magnitude greater than those of 
the clusters. 



The structure of both objects is similar to that of IRAS 
23385+6053. This source, which has a bolometric lumi- 
nosity of ^ 10** Lq, has been proposed by Molinari et al. 
( 1998bj and Fontani et al. ( 2004) as the best example of a 
candidate massive protostar. As in 121307 and 122172, we 
can distinguish two components: a compact core detected 
at millimeter and sub-millimeter wavelengths and a clus- 
ter of stars which lie around the core, detected at near- 
and mid- infrared wavelengths. Other important similari- 
ties between IRAS 23385+6053, 121307 and 122172 are the 
lack of free-free emission towards the core and the pres- 
ence of a compact outflow driven by a source embedded 
in the core. 

Given the strong similarities between 121307, 122172 
and IRAS 23385+6053, it seems plausible that also the 
nature of the embedded YSOs is approximately the same. 
Molinari et al. Ijl998b|l propose two possible scenarios for 
IRAS 23385+6053, that the object which lies inside the 
core is a BO ZAMS star, and the Hii region is not de- 
tectable because it is very compact and optically thick; 
another possibility is that there is no Hii region because 
the central object is still in the proto-stellar phase. In 
both scenarios, the embedded source is young and mas- 
sive, while the stars of the surrounding cluster are more 
evolved. These scenarios are plausible also for 121307 and 
122172. In the first hypothesis, the objects embedded in- 
side I21307-C and I22172-C are B1-B2 ZAMS stars; in the 
second, they are intermediate-mass protostars, and the lu- 
minosity observed is due to accretion. Behrend & Maeder 
(I 2001II computed evolutionary tracks for protostars from 
1 to 85 Mq: in their models the luminosities of our sources 
correspond to protostellar masses of ~ 5 — 8 Mq. However, 
this value must be regarded with caution, because the au- 
thors assume the relation found by Churchwell (1998) be- 
tween L and Mout, which in the high-mass regime does 
not provide a good fit to the data. 

As discussed in Sect. 14.21 the protostellar hypothesis 
seems to be more likely. However, to confirm this conclu- 
sion further interferometric observations are needed. In 
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particular, it would be important to measure the tem- 
perature of the cores with great accuracy. In fact, newly 
formed massive stars lie inside "hot cores", namely com- 
pact molecular regions with kinetic temperature higher 
than ~ 100 K (Kurtz et al. 2000) . Hence, as discussed 
by Fontani et al. I)2U(J4|I . measuring the core temperature 
would be an important step to improve our understanding 
of their evolutionary stage. 

4.3. The formation of "Cluster! nos" 

As we have discussed in the previous sections, it is very 
likely that the outflows, the bulk of the bolometric lumi- 
nosity and other properties of the regions are dominated 
by the most massive forming star. Nevertheless, our in- 
frared images reveal that small young stellar aggregates 
are also found in the same regions (see also Molinari et 
al. 1998b and Fontani et al. 2004). It is interesting to ex- 
plore the idea that, along with the 5-8 Mq protostar, 
a small stellar cluster is forming within the molecular 
clump, a "clusterino" similar to those found around the 
more evolved early type Herbig Be stars (see e.g. Testi et 
al. 11999(1 . The values of the clusterinos richness indicators 
that we derived in Sect. 13. 2 H are relatively low if compared 
with the values for Herbig Be clusters. This is not unex- 
pected, as our observations are not as sensitive (in terms 
of the least massive objects detectable) as those of Testi 
et al. H1998fl . this is mainly due to the larger distance of 
our sources compared with the sample of Herbig Be stars. 
An additional reason could be that our objects are gener- 
ally younger and more extincted, and, finally, it is possible 
that the clusterinos have not been completely assembled 
yet. 

From the photometric data alone, it is difficult to es- 
timate the membership of individual objects (see also the 
discussion in the Testi et al. papers on Herbig Ae/Be clus- 
terinos). Assuming that all the sources within ~30" from 
the IRAS nominal position are cluster members and all 
the other sources in the field are not related to the clus- 
terinos, it is possible to derive an estimate of the most 
luminous (already formed) star and the integrated near in- 
frared emission for each clusterino. The result of this exer- 
cise shows that the most massive star in the 121307 region 
is consistent with a moderately reddened (Ay ~ 4mag) 
B6 ZAMS, while in 122172 the brightest near infrared star 
close to the IRAS position could be a more embedded 
(Ay ~ 9mag) B1.5 ZAMS. The integrated near infrared 
fluxes at ~2.2 /im are ^5 mJy and ^^90 mJy for 121307 
and 122172, respectively. These numbers show that the al- 
ready formed members of the clusterinos contribute only 
to a minimal fraction of the emission from the region in 
the case of 121307, while in 122172 a substantial fraction of 
the total luminosity could be ascribed to the most massive 
member. 

It is important to point out that these estimates of 
the most massive member of the clusterinos and their in- 
tegrated emission at near infrared wavelengths have been 



made assuming that all objects close to the IRAS source 
are bona fide members. Clearly this is a simplifying as- 
sumption: it is likely that an important contribution from 
foreground and background objects is affecting our es- 
timates. The comparison of the richness indicator val- 
ues and the number of objects close to the IRAS source 
shows that this contamination could be as high as 50%. 
Unfortunately, with our data alone, it is not possible to 
pin point the contaminating objects, hence our estimates 
should be regarded as upper limits (both for the luminos- 
ity of the most massive star and the integrated flux of the 
clusterinos). 

5. Summary and conclusions 

We have observed two massive proto-stellar candidates, 
IRAS 21307+5049 and IRAS 22172+5549, in the ^^CO (1- 
0) and (2-1), C^^O (2-1), Hi3C0+ (1-0), CH3C2H (6-5), 
(8-7) and (13-12) molecular lines, and in the continuum 
at various wavelengths from the near-infrared to the mil- 
limetric range. The following results have been obtained: 

— Maps of the 3 mm continuum observed with the OVRO 
interferometer reveal in both sources the presence of a 
compact, massive core, whose linear diameter and mass 
are - 0.09 pc and - 50 Mq for 121307, and - 0.04 pc 
and - 40 Mq for 122172; 

— Near- infrared images at 1.7 and 2.2 /im indicate that 
clusters of more evolved stars surround both cores. The 
continuum maps and the SEDs suggest that for both 
sources the core is mainly responsible for the emission 
at wavelengths longer than ~ 25 /im, while the emis- 
sion at shorter wavelengths is due to the surrounding 
cluster. The SEDs indicate that the core luminosity is 
a few 10'^ Lq for both sources. 

— By mapping the ^^CO (1-0) and (2-1) line wings, 
we find outflows associated with both sources that 
are driven by YSOs embedded in the cores. Outflows 
have masses of a few to a few tens of Mq, dynami- 
cal timescales of ^ W^ — 10^ yr and mass loss rates 
of '^ a few times ^ lO~^M0yr~^. Images of the H2 
and [Fell] lines at near-infrared wavelengths suggest 
the presence of an ionized jet in the inner part of the 
outflow associated with 122172. 

— The flows physical parameters of both sources are con- 
sistent with those of a sample of massive protostar 
candidates (Beuther et al. 2pQ2a) with comparable lu- 
minosity. The values of the accretion efficiencies, of 
~ 10~^, are consistent with values derived assum- 
ing that in these sources the luminosity is accretion- 
dominated. 

— Based on the similarity between our sources and the 
massive protostar candidate IRAS 23395+6053, we be- 
lieve that the massive YSOs embedded inside the cores 
are protostars accreting material from the environ- 
ment, the masses of which at present are ~ 5 — 8 Mq, 
assuming that the luminosity observed is due to accre- 
tion onto the forming star. 
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— Near-infrared images suggest that, together with the 
massive protostar, a small stellar cluster ( "clusterino" ) 
is forming inside the molecular clump, similar to those 
found around Herbig Be stars. Assuming that the 
sources within 30" from the IRAS position are all bona 
fide members of the clusterino, we have derived that 
the already formed stars of the clusterino contribute 
only to a small fraction to the total observed luminos- 
ity in 121307, while in 122172 a significative fraction of 
the total luminosity could be due to the most massive 
members. 
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Fig. 14. Same as Fig.^lfor 122172. a: contour levels of the 
blue- and red- wings begin from 8 ('~4a) K km s""'^, with a 
spacing of 4 K km s^^. The velocity ranges are (—51,-45) 
km s~^ (blue) and (-38,-30) km s~^ (red). The grey- 
scale represents the SCUBA 850 /im emission map: the 
same contour levels as in Fig. |S1 are shown, b: contour 
levels begin from 5 K km s"^ (~ 3tT level in the maps) and 
have a spacing of 5 K km s~^. c: same as c. of Fig. 1131 
Contour levels of the blue and red wings begin from 1.2 
K km s~^, and have a spacing of 0.4 K km s~^. The grey- 
scale represents the OVRO 3 mm map, identifying the 
core I22172-C, with the same contours as in Fig.|Sl 
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Fig. 16. Mechanical force of the outflows, Poutj versus 
bolometric luminosity Lboi of the source. 121307 and 
122172 data are represented by squares and pentagons, 
respectively (open symbol = OVRO data, filled symbols 
— single-dish data); crosses indicates data from Bontemps 
et al. (|lH9(i(l (low-mass YSOs) and open circles correpond 
to those of Beuther et al. (2002a) (high-mass YSOs). The 
filled triangle indicates IRAS 23385+6053 (Fontani et al. 
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Fig. 15. [Fell] and H2 line images of 122172 after contin- 
uum subtraction. The solid and dashed lines indicated the 
blue and red lobes of the outflow, respectively, traced by 
the ^^CO (1-0) line. The same contour levels as in Fig. 1141 
are shown. The position of the core I22172-C is indicated 
by the thick contour, which represents the half maximum 
power level of the OVRO map. 
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Fig. 17. Accretion efficiencies estimated from Eq. (O 
(crosses), from Eq. (jSJ) (filled circles) and from Eq. JTJ 
(open circles) , against the bolometric luminosity. Big sym- 
bols correspond to 121307, 122172 and IRAS 23385+6053; 
small symbols indicate the sources of the Beuther et al. 
I|2002a|l sample. 



